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Abstract 

A  finite  element  method-based  procedure  is  developed  for  the  design  of  molecularly  functionalized  nano-size  devices.  The  procedure  is 
aimed  at  the  single- walled  carbon  nano-tubes  (SWCNTs)  used  in  the  construction  of  such  nano-devices  and  utilizes  spatially  varying  nodal 
forces  to  represent  electrostatic  interactions  between  the  charged  groups  of  the  functionalizing  molecules.  The  procedure  is  next  applied  to 
the  design  of  a  fluid-flow  control  nano-valve.  The  results  obtained  suggest  that  the  finite  element-based  procedure  yields  the  results,  which  are 
very  similar  to  their  molecular  modeling  counterparts  for  small-size  nano- valves,  for  which  both  types  of  analyses  are  feasible.  The  procedure 
is  finally  applied  to  optimize  the  design  of  a  larger-size  nano- valve,  for  which  the  molecular  modeling  approach  is  not  practical. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  the  last  several  years,  considerable  effort  has  been  in¬ 
vested  into  studying  functionalized  (typically  silicon-based) 
cantilever  nano-beams  e.g.  [1,2].  Such  nano-beams  have  one 
face  covered  (functionalized)  with  self-assembled  or  chem¬ 
ically  bonded  molecules  of  various  types  and  can  deflect  by 
themselves  due  to  differences  in  the  surface  stress  between 
the  functionalized  and  non-functionalized  faces.  In  a  recent 
paper,  Solares  et  al.  [3  ]  proposed  the  design  of  a  fluid-flow 
control  nano-valve  based  on  a  silicon  nano-beam  actuator 
functionalized  with  a  covalently  bonded  monolayer  of  acrylic 
acid  molecules,  which  can  allow  or  stop  the  flow  through  a 
fluid  conduit,  a  single-walled  carbon  nano-tube  (SWCNT). 
The  “on/off’  switching  of  the  nano-valve  is  controlled  by  the 
pH  of  its  environment.  That  is,  changes  in  the  pH  affect  the 
charges  of  the  bonded  molecules  and,  thus,  control  the  mag¬ 
nitude  of  the  inter-molecular  electrostatic  forces. 

Solares  et  al.  [3]  argued  that  as  molecular  modeling  is 
moving  from  being  primarily  a  tool  for  studying  materials 
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properties  to  an  alternative  approach  for  analyzing  the  per¬ 
formance  of  devices,  it  is  becoming  computationally  more 
and  more  demanding  and,  often,  impractical.  This  is  caused 
by  the  fact  that  the  size  of  the  computational  problem  gen¬ 
erally  increases  as  N2,  where  N  is  the  number  of  particles  in 
the  system  being  studied.  This  problem  is  particularly  aggra¬ 
vated  by  the  presence  of  electrostatic  interactions  between 
the  charged  molecules  used  for  nano-beam  functionalization 
because  of  the  following  two  main  reasons:  (a)  the  geome¬ 
try  of  nano-devices  is  generally  not  periodic  so  that  efficient 
electrostatic  lattice-summation  methods  (e.g.  the  Ewald  sum¬ 
mation  method)  cannot  be  used  without  introducing  spurious 
effects  due  to  imposition  of  the  periodic  boundary  conditions 
and  (b)  the  use  of  radial  cut-off  radii  or  splines  to  handle 
the  long-range  interactions  can  underestimate  the  electro¬ 
static  energy  by  as  much  as  an  order  of  magnitude  for  the 
nanometer-size  devices.  There  is  an  additional  problem  asso¬ 
ciated  with  the  use  of  molecular  modeling  methods  for  sim¬ 
ulation  of  the  nano-size  devices.  Average  forces  and  strain 
energies  used  to  judge  the  attainment  of  equilibrium  may 
not  be  very  reliable  in  a  computational  system  containing  a 
large  number  of  particles,  since  they  may  allow  large  resid¬ 
ual  forces  to  be  present  in  a  small  region,  while  the  average 
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force  is  quite  small  implying  an  equilibrium  state  of  the  sys¬ 
tem. 

Nevertheless,  to  obtain  a  more  accurate  design  of  a  nano¬ 
device,  it  is  critical  to  include  as  much  molecular  level  details 
as  possible.  As  pointed  out  above,  on  the  other  hand,  it  is  often 
impractical  or  impossible  to  carry  out  all  design  and  analy¬ 
sis  steps  at  the  molecular  level.  To  achieve  a  compromise, 
Solares  et  al.  [3]  proposed  a  hybrid  method,  which  com¬ 
bines  molecular  simulations  with  the  classical  engineering 
design.  Specifically,  molecular  simulations  are  used  to  de¬ 
termine  elastic  and  electrostatic  properties  of  the  individual 
components  of  a  nano-device,  which  are  then  used,  within 
a  classical  engineering  design  framework,  to  carry  out  an 
analysis  of  the  functionality,  design,  and  optimization  of  the 
nano-device. 

The  objective  of  the  present  work  is  to  extend  the  hybrid 
method  proposed  by  Solares  et  al.  [3]  and  apply  it  to  the  de¬ 
sign  of  a  fluid-control  valve  described  above.  Specifically,  the 
method  of  Solares  et  al.  [3]  fails  to  include  the  interactions 
between  different  device  components  (primarily  the  interac¬ 
tions  between  the  cantilever  nano-beam  and  the  SWCNT) 
during  the  operation  of  such  a  nano- valve.  Rather,  the  perfor¬ 
mances  of  the  cantilever  nano-beam  and  the  SWCNT  fluid 
conduit  are  analyzed  separately  and  that  information  used  to 
judge  the  functionality  of  the  nano-device  and  determine  its 


optimal  geometry  and  dimensions.  In  the  present  work,  we 
extend  the  hybrid  method  proposed  by  Solares  et  al.  [3]  by 
utilizing  a  finite  element  method  to  analyze  the  performance 
of  the  entire  nano-valve  and  use  the  results  of  such  analysis 
as  a  guide  for  optimization  of  the  device  design.  In  this  way, 
the  interactions  between  the  cantilever  nano-beam  and  the 
SWCNT,  which  are  at  the  heart  of  operation  of  the  function¬ 
alized  silicon  cantilever  nano-beam-based  nano-valve,  are  di¬ 
rectly  taken  into  account. 

The  organization  of  the  paper  is  as  follows:  a  brief  for¬ 
mulation  to  the  problem  being  studied  is  given  in  Section 
2.1.  The  use  of  the  shell  theory  to  model  the  SWCNT  wall 
within  a  continuum  formulation  is  presented  in  Section  2.2. 
The  development  of  special  contact  elements,  which  account 
for  the  van  der  Waals  interactions  between  contacting  sur¬ 
faces  is  presented  in  Section  2.3.  The  procedure  used  to  com¬ 
pute  energy  and  force  associated  with  the  electrostatic  in¬ 
teraction  between  the  charged  functionalizing  molecules  is 
explained  in  Section  2.4.  The  method  employed  to  specify 
the  initial  internal  stress  in  the  SWCNT  wall  is  described  in 
Section  2.5.  A  brief  overview  of  the  finite  element  method 
used  is  given  in  Section  2.6.  The  main  results  obtained  in  the 
present  work  are  presented  and  discussed  in  Section  3,  while 
the  key  resulting  conclusions  are  summarized  in  Section 
4. 


SWCNT  Fluid  Conduit 


Monolayer  of  Functionalizing  Molecules 


Fig.  1.  Atomistic  structure  of  the  fluid-flow  control  nano- valve  analyzed  in  the  present  work  in  the:  (a)  open  and  (b)  closed  positions. 
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2.  Computational  procedure 

2.1.  Problem  definition 

Atomistic  configurations  of  the  nano-valve  in  the  open 
and  closed  positions  are  displayed  in  Fig.  1(a)  and  (b),  re¬ 
spectively.  The  nano-valve  consists  of  a  silicon  block  with  a 
cantilever  nano-beam,  where  both  are  functionalized  on  their 
top  surface  with  a  single-layer  of  molecules.  The  silicon  block 
is  perforated  on  its  lower  part  and  a  SWCNT  inserted  through 
the  perforation.  The  SWCNT  acts  as  a  fluid  conduit.  When  the 
pH  in  the  environment  surrounding  the  nano-valve  changes, 
the  functionalizing  molecules  may  acquire  a  larger  charge. 
This,  in  turn,  gives  rise  to  a  larger  electrostatic  repulsion  be¬ 
tween  the  molecules  and,  subsequently,  causes  the  cantilever 
to  bend  downward.  If  the  action  of  the  cantilever  causes  the 
SWCNT  to  deflect  past  the  point  of  bucking,  the  fluid-flow 
through  the  SWCNT  ceased  and  the  nano-valve  is  in  the  "off' 
position.  Fig.  1(b). 

For  the  nano-valve  to  function  as  described  above,  the  re¬ 
duction  in  the  electrostatic  energy  associated  with  an  increase 
in  the  inter-molecular  separation  during  cantilever  bending 
has  to  be  sufficiently  large  to  compensate  for  the  elastic  ener¬ 
gies  stored  in  the  cantilever  nano-beam  and  in  the  nano-tube 
(up  to  and  past  the  point  of  buckling  of  the  SWCNT).  For 
a  given  charge-level  of  the  functionalizing  molecules,  the 
electrostatic  energy  reduction  is  governed  by  the  molecular 
surface  coverage  (density)  and  the  surface  area  of  the  top  face 
of  the  block  and  the  cantilever.  For  a  fixed  molecular  surface 
density  and  a  fixed  value  of  the  combined  lengths  of  the  sili¬ 
con  block  and  the  cantilever,  the  nano-valve  width  becomes 
the  sole  design  parameter  controlling  the  functionality  of  the 
nano-valve. 

2.2.  Application  of  shell  theory  to  SWCNTs 


wise,  the  bending  stiffness  is  mainly  a  measure  of  the  resis¬ 
tance  to  C— C— C  (two-bond)  in-plane  bond-angle,  C—C—C—C 
C 

dihedral  torsion-angle  and^,  >  C  -  C  out-of-plane  inversion- 
angle  changes. 

Following  the  recommendation  of  Pantano  et  al. 
[4],  representative  values  for  the  stretching  stiffness, 
C  =  59.36eV/atom  =  400.4  J/m2,  and  for  the  bending  stiff¬ 
ness,  D  =  2.886  eV  A2/atom=  1.765  x  10“ 19  J  are  used.  The 
effective  SWCNT  wall  Young’s  modulus,  £wain  and  the  ef¬ 
fective  SWCNT  wall  thickness,  ?wan,  are  computed  from  the 
following  relations: 

C  =  £walPwall^2  (1) 

and 

D  =  (2) 

12d  -  ''wall) 

as  Fwaii =  4.84  TPa  and  fwaii  =  0.075  nm,  where 
£2  =  0.0262  nm2/atom  is  the  surface  area  of  the  SWCNT  wall 
allocated  to  a  carbon  atom,  and  i>waii  =  0. 1 9  is  the  SWCNT 
wall  Poisson’s  ratio  (set  equal  to  that  in  a  graphene  sheet). 

The  SWCNT  wall  (i.e.  the  cylindrical  shell  wall)  is  dis¬ 
cretized  within  the  finite  element  framework  used  in  the 
present  work  using  four-node  doubly  curved  shell  elements, 
which  allow  transverse  shear  deformation.  These  elements 
are  suitable  for  modeling  the  SWCNT  wall  using  a  contin¬ 
uum  elastic  shell  model,  which  accounts  for  geometrical  non- 
linearities  (critical  for  modeling  the  onset  of  buckling  in  the 
SWCNT).  The  reference  surface  of  a  (SWCNT  wall)  shell 
is  defined  in  terms  of  the  shell  element  nodal  coordinates 
and  the  normal  definitions.  The  shell’s  degrees  of  freedom 
including  all  kinematic  quantities  are  all  referred  to  such  a 
reference  surface. 


Following  previous  work  e.g.  [4],  a  SWCNT  is  treated 
within  a  continuum  approach  as  a  cylindrical  thin-wall  shell. 
The  stretching  and  the  bending  stiffnesses  of  such  a  shell 
can  be  determined  using  ab  initio  or  atomistic  simulations 
of  the  stretching  and  the  bending  behavior  of  SWCNTs  e.g. 
[4].  Using  standard  relations  given  below,  the  two  stiffnesses 
can  next  be  converted  into  the  effective  elastic  moduli  (the 
Young’s  modulus  and  the  Poission’s  ratio)  and  the  effective 
thickness  of  the  SWCNT  (shell  wall). 

Within  an  atomistic  simulation  framework  within  which 
bonding  between  the  carbon  atoms  in  a  SWCNT  is  repre¬ 
sented  using  a  forcefield,  all  terms  in  the  forcefield  potential 
affect  the  resistance  of  a  SWCNT  to  any  mode  of  deforma¬ 
tion.  However,  the  stretching  stiffness  is  generally  dominated 
by  the  resistance  of  SWCNTs  to  C— C  bond  stretching  and  to 
C—C—C  (two-bond)  in-plane  bond-angle  changes.  Within  the 
condensed-phased  optimized  molecular  potential  for  atom¬ 
istic  simulation  studies  (COMPASS)  forcefield  [11],  for  ex¬ 
ample,  these  two  forcefield  terms  are  responsible  for  over 
95%  of  the  contribution  to  the  stretching  resistance.  Like- 


2.3.  Development  of  special  interaction  elements 


To  further  include  atomic-level  effects  into  the  present 
continuum-level  analysis,  van  der  Waals  interactions  between 
neighboring  surfaces  must  be  taken  into  account.  Specifically, 
the  van  der  Waals  interactions  between  different  sections  of 
the  internal  surface  of  a  SWCNT  must  be  considered,  since 
such  interactions  affect  the  post-buckling  behavior  and  the 
reversibility  of  large-scale  deformation  in  a  SWCNT.  In  ad¬ 
dition,  the  van  der  Waals  interactions  between  the  tip  of  the 
silicon  cantilever  nano-beam  and  the  SWCNT  fluid  conduit 
control  the  transfer  of  the  force  from  the  nano-beam  to  the 
SWCNT,  and,  hence,  must  be  taken  into  account. 

The  van  der  Waals  interactions  within  a  SWCNT  are 
modeled  using  the  following  Lennard-Jones-based  pres¬ 
sure/interlayer  distance  relation  [4]: 


Pc-c 


I  (!) 


(3) 
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where  Pq~c  is  the  carbon-carbon  interlayer  contact  pres¬ 
sure,  a  the  interlayer  distance,  do  =  0.34  nm,  the  equilibrium 
graphene  inter  plane  distance  and  tfr  =  36.5  GPa. 

To  obtain  a  similarpressure/interlayer  distance  relation  for 
the  nano-beam/SWCNT  interactions,  a  short  molecular  me¬ 
chanics  simulation  analysis  is  carried  out  using  the  Accelrys’ 
molecular  modeling  program.  Discover  [5],  and  the  COM¬ 
PASS  forcefield  potential  [6].  Since  the  silicon  block  and  the 
cantilever  beam  used  in  the  present  work  have  a  cubic  (001) 
orientation,  a  computational  domain  is  constructed,  which 
consists  of  one  (00  1)  silicon-plane  and  one  graphene  plane. 
The  in-plane  size  of  the  computational  domain  in  terms  of 
the  number  of  silicon  lattice  parameters  used  is  seven  by 
nine.  The  size  of  the  corresponding  graphene  sheet  was  dif¬ 
ferent  by  0.9  and  0.7%  in  the  x-  and  y-direction  from  those  of 
the  silicon-plane  (00  1).  To  eliminate  this  difference,  equal- 
magnitude/opposite-sign  strains  are  applied,  respectively,  to 
the  silicon  and  the  graphene  sheets,  in  both  directions.  The 
periodic  boundary  conditions  are  next  applied  in  the  x-  and 
y-direction.  The  initial  inter  planar  distance  was  then  set  to 
0.8  nm  and  reduced  repeatedly  in  decrements  of  0.01  nm.  At 
each  value  of  the  inter  planar  spacing,  the  interaction  energy 
of  the  two  planes  is  computed.  The  interaction  energy  is  de¬ 
fined  as  a  difference  between  a  sum  of  the  potential  energies 
of  the  two  planes  at  an  infinite  distance  and  the  energy  of 
the  computational  domain  containing  the  two  planes  at  a  fi¬ 
nite  distance.  The  results  of  this  calculation,  as  provided  in 
Fig.  2(a),  are  fitted  to  the  following  functional  form: 


Pint  —  C 


■nfdoy 

jn  \  a  J 


(4) 


where  m  and  n  <  m  are  positive  integers,  and  C  is  a  positive 
proportionality  constant. 

The  Levenberg-Marquart  method  for  non-linear  parame¬ 
ter  estimation  is  next  used  to  determine  the  unknown  param¬ 
eters:  m,  n,  C  and  do  in  Eq.  (4).  The  fitting  procedure  yielded 
the  following  optimal  values  for  the  model  parameters:  m  =  l , 
n  =  4,  C  =  0.6539  J/m2  and  do  =  0.364  nm.  The  results  of  the 
fitting  procedure  are  also  displayed  in  Fig.  2(a).  It  is  seen 
that  Eq.  (4)  can  quite  successfully  account  for  the  interaction 
energy/inter  planar  spacing  relation  obtained  using  atomistic 
simulations. 

Eq.  (4)  is  next  used  along  with  the  values  of  the  optimal 
parameters  to  define  the  pressure/interlayer  distance  as: 


Fsi-c  =  - 


dgint 

d  a 


Cn  T /do \,,,+1 
do  V  «  / 


(5) 


where  Psi-C  in  the  Si(0  0  l)/graphene  interlayer  contact  pres¬ 
sure. 

The  pressure/interlayer-spacing  relations  given  by  Eqs.  (3) 
and  (5)  are  implemented  in  the  ABAQUS  user  subroutine 
UINTER  [7],  which  enables  the  user  to  specify  the  constitu¬ 
tive  interactions  between  two  deforming  surfaces  (a  master 
and  a  slave  surface).  The  subroutine  is  called  for  each  slave 
node  at  each  time  increment  during  the  finite  element  anal¬ 


Inter-planar  Distance,  nm 

Fig.  2.  Variation  of:  (a)  the  Si(0  0  l)/graphene  interaction  energy  and  (b)  the 
contact  pressure  with  the  inter  planar  separation. 


ysis  to  compute  the  tractions  at  such  nodes  based  on  the  in¬ 
crements  in  their  relative  position  with  respect  to  the  master 
surface.  The  normal  tractions  can  be  either  positive  (indicates 
surfaces  repulsion)  or  negative  (denotes  surfaces  attractive). 
In  addition  to  computing  the  nodal  traction,  the  UINTER  sub¬ 
routine  is  used  to  compute  the  corresponding  Jacobian  to  help 
to  accelerate  convergence  of  the  computational  algorithm. 

2.4.  Evaluation  of  the  monolayer  electrostatic  energy 
and  force 

The  electrostatic  energy  and  force  for  a  functionalizing 
single-layer  can  be  evaluated  by  placing  charged  particles  at 
the  positions  on  the  top  face  of  the  silicon  block  and  can¬ 
tilever  nano-beam,  where  the  charged  functional  groups  of 
the  functionalizing  monolayer  would  be  located  during  bend¬ 
ing.  Such  a  procedure  yields  a  functional  form  between  the 
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Nano-beam,  nm 

Fig.  3.  Variations  of  the  initial  nodal  forces  for  different  nano-valve  widths 
(dashed  lines)  and  the  initial  minimal  requires  nodal  forces  (solid  line)  along 
the  length  of  the  cantilever  nano-beam  (measured  from  the  fixed  end  of  the 
nano-beam). 

total  electrostatic  energy  or  the  electrostatic  forces  acting  on 
each  charged  particles  at  a  given  level  of  bending  curvature 
of  the  cantilever.  Ideally,  one  would  like  to  use  such  forces 
within  a  static  finite  element  analysis  like  the  one  used  in 
the  present  work  (described  in  the  next  section).  However, 
the  finite  element  program  ABAQUS  [7]  used  in  the  present 
work  does  not  allow  the  applied  force  to  be  varied  during  the 
simulation  run,  i.e.  only  time-independent  concentrated  or 
distributed  forces  can  be  used.  To  overcome  this  limitation, 
the  procedure  described  below  was  utilized,  which  involves 
the  following  steps: 

(a)  A  cantilever  nano-beam  with  a  sufficient  width  is  first 
selected  to  ensure  that  the  electrostatic  forces  acting  be¬ 
tween  the  charged  particles  are  sufficiently  large  to  give 
rise  to  the  nano-tube  buckling.  For  a  give  molecular  sur¬ 
face  coverage  and  particle  charges,  the  force  on  each 
charged  particle,  i,  due  to  its  interaction  with  all  other 
particles,  j,  is  next  computed  using  the  following  equa¬ 
tion: 


(b)  Once  such  forces  are  computed  for  all  the  nodes  at  the 
top  face  of  the  cantilever  nano-beam,  a  finite  element 
analysis  is  carried  out  until  the  onset  of  SWCNT  buck¬ 
ling.  At  that  point,  the  coordinates  of  all  the  nodes  at  the 
top  face  of  the  cantilever  nano-beam  are  determined  and 
the  forces  acting  on  each  corresponding  charged  particle 
calculated  using  Eq.  (6).  Since  the  charged  particles  are 
separated  further  apart  in  the  deformed  cantilever  con¬ 
figuration,  the  forces  computed  are  lower  in  magnitude 
relative  to  their  counterparts  in  the  original  nano-tube 
configuration. 

(c)  Next,  it  is  imposed  that  the  work  done  by  each  node 
under  the  condition  of  a  constant  force  between  the  ini¬ 
tial  and  the  final  configuration  (as  encountered  during 
a  finite  element  method  run)  is  identical  to  that  which 
would  be  done,  if  the  force  acting  on  that  node  were 
allowed  to  vary  with  the  nano-beam  curvature.  This  re¬ 
quirement  and  the  fact  that  the  electrostatic  force  de¬ 
creases  as  1  /rfj  during  cantilever  bending  were  then 
used  to  determine  the  appropriate  initial  force,  which 
should  act  on  each  of  the  nodes.  A  variation  of  such 
a  minimum  required  initial  force  acting  on  the  nodes 
with  a  distance  along  the  cantilever  nano-beam  for  a 
nano-valve  with  a  silicon  block  length  of  40,  the  can¬ 
tilever  nano-beam  length  of  30  and  the  width  of  14  (all 
in  terms  of  the  number  of  silicon  lattice  parameters)  is 
shown  in  Fig.  3,  curve  labeled  “minimal  required  initial 
force”. 

(d)  The  variation  of  the  force  acting  on  the  individual  nodes 
along  the  nano-beam  length  in  the  initial  configuration 
is  next  generated  using  the  procedure  described  in  part 
(a)  for  several  nano-beam  widths.  The  nano-beam  widths 
are  incremented  by  an  integer  number  of  silicon  lattice 
parameters.  The  smallest  nano-beam  width,  which  sat¬ 
isfies  the  condition  that  the  force  acting  at  any  node  in 
the  initial  configuration  is  greater  than  the  minimal  re¬ 
quired  initial  force  is  deemed  as  the  minimal  acceptable 
cantilever  beam  width. 


2.5.  Initial  internal  state  within  the  SWCNT 


F  i  = 


■y-'  Ritfj . 


4jT£o 


(j  #  0 


(6) 


where  q,  and  qj  are  the  charge  of  particles  i  and  j,  rq  the  i-j 
inter  distance,  rq  the  unit  vector  point  along  the  direction 
connecting  particle  i  to  j,  so  the  permittivity  of  free  space 
and  bold  fonts  are  used  to  denote  vector  quantities. 
Since  the  plane-strain  boundary  conditions  are  applied 
to  the  cantilever  nano-beam  in  the  width  direction,  the 
forces  acting  on  all  atoms  at  a  given  location  along  the 
cantilever  length  (i.e.  at  a  given  .r-coordinate)  are  first 
summed,  then  divided  by  2  and  finally  assign  to  each  of 
the  two  nodes  at  a  given  x-location  at  the  top  face  of  the 
cantilever  nano-beam  Fig.  3. 


As  pointed  out  by  several  research  groups  e.g.  [8,9],  the 
SWCNTs  contain  an  internal  stress  due  to  their  curvature. 
Pantano  et  al.  [4]  observed  that  such  internal  stress  can  play 
a  major  role  in  the  ability  of  a  SWCNT  to  recover  large  lo¬ 
calized  strains  and,  hence,  must  be  taken  into  account  in  the 
present  case,  where  such  strain  reversibility  is  critical  aspect 
of  the  “ on/off ’  switching  behavior  of  the  fluid-flow  nano¬ 
valve.  To  compute  the  initial  circumferential  internal  stress 
within  the  SWCNT,  erg,  the  rolling  (strain)  energy  of  a  single 
graphene  sheet: 


Cr  = 


D 


9  e>2 
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is  differentiated  with  respect  to  nano-tube  radius,  /?swcnt> 
to  get: 

a9  =  -  'wan^wan - 1  (8) 

12^  ''wall)  ^SWCNT 

where  /cWaii  =  1  /Rswcnt  is  the  SWCNT  wall  curvature. 

The  circumferential  internal  stress  given  by  Eq.  (8)  is  in¬ 
corporated  into  the  ABAQUS  subroutine  SIGINI  in  order  to 
specify  the  initial  stress  state  in  the  nano-tube  (shell)  wall. 

2.6.  Finite  element  analysis 

While  many  commercial  finite  element  codes  are  capa¬ 
ble  of  handling  geometrically  non-linear  small-strain/large- 
rotation  shell  (used  to  represent  the  SWCNTs)  and  brick 
(used  to  represent  silicon  cantilever  nano-beam)  elements, 
ABAQUS  standard  [7]  is  chosen  in  the  present  investigation, 
as  indicated  earlier,  since  it  contains  user  subroutines  for  im¬ 
plementation  of  the  van  der  Waals  interactions  and  the  initial 
internal  stress  state  in  the  nano-tube. 

A  schematic  representation  of  the  finite  element  model  is 
shown  in  Fig.  4  along  with  the  appropriate  initial  and  bound¬ 
ary  (loading)  condition.  The  SWCNT  is  represented  using 
four-node  doubly  curved  reduced-integration  elements  S4R 
with  hourglass  control.  The  silicon  cantilever  is  represented 
using  eight-node  brick  element  C3D8  ABAQUS  elements. 
Following  Pantano  et  al.  [4],  the  S4R  shell  elements  are  se¬ 
lected  to  be  approximately  square  in  shape  and  to  have  a  side 
length  around  0.24  nm,  the  height  of  the  hexagonal  cell  in 
the  graphene  plane.  The  size  of  the  C3D8  elements  in  the  x- 
and  v-direction  is  chosen  to  nearly  match  that  of  the  shell  ele¬ 
ments.  Almost  identical  results  of  the  finite  element  analyses 
are  obtained  using  elements,  which  are  larger  or  smaller  by 
a  factor  of ,  suggesting  that  a  mesh-size  invariance  condition 
has  been  attained. 

The  SWCNT  material  is  modeled  as  an  elastically 
isotropic  material  with  Fwa\\  and  vwau  given  in  Section  2.2. 
The  silicon  material  is  also  modeled  as  an  isotropic  lin¬ 
early  elastic  material  with  Fbeam  =  76.7  GPa  and  Vbeam  =  0.26. 
These  values  were  obtained  by  fitting  the  bending  en¬ 
ergy  results  obtained  using  atomic-level  bending  simula¬ 
tions  to  a  standard  bending  energy /bending  curvature  relation 

[3]. 


Fig.  4.  The  formulation  of  the  finite  element  boundary  value  problem  ana¬ 
lyzed  in  the  present  work. 


3.  Results  and  discussion 

3.1.  Validation  of  the  present  finite  element  method 

For  the  silicon  cantilever  beam-based  fluid-flow  nano¬ 
valve  analyzed  in  the  present  work  to  be  functional,  the  di¬ 
ameter  of  the  SWCNT  fluid  conduit  must  be  large  enough 
to  enable  a  free  passage  of  the  molecules,  which  as  in  the 
case  of  organic  substances  or  drugs,  may  have  a  large  size. 
A  large  diameter  SWCNT,  however,  requires  a  large  sili¬ 
con  block,  a  large  cantilever  nano-beam  and,  thus,  a  large 
number  of  functionalizing  molecules  in  the  top  face  mono- 
layer.  Consequently,  the  number  of  particles  in  the  device 
may  become  prohibitively  large  for  the  problem  to  be  ana¬ 
lyzed  using  molecular  modeling.  That  is,  the  reason  why  a 
finite  element  method-based  approach,  which  does  not  suffer 
from  the  similar  limitations,  is  developed  in  the  present  work. 
However,  the  molecular  modeling  approach  is  intrinsically 
more  accurate  and  reliable  than  its  finite  element  counterpart 
and,  hence,  one  must  develop  a  procedure  for  validating  the 
latter. 

The  procedure  used  in  the  present  work  to  validate  the 
finite  element-based  approach  involves  finding  the  solution 
to  the  same  problem  by  using  both  a  more  accurate  molec¬ 
ular  modeling  method  and  the  less  accurate  finite  element 
method-based  approach  presented  in  Section  2.6.  Due  to  the 
limitations  of  the  molecular  modeling  approach  described 
above,  a  nano-valve  of  a  smaller  size  than  the  one,  which 
could  function  as  a  regulator  for  the  flow  of  organic/drug 
molecules  are  analyzed.  Specifically,  a  (8,  8)  SWCNT  (the 
nano-tube  radius  =  0.54  nm,  the  nano-tube  length  =  12.0  nm) 
is  selected  as  the  flow  conduit  and  the  lengths  of  the  silicon 
block  and  the  cantilever  length  are  both  set  to  5  nm. 

Next,  the  procedure  described  in  Sections  2.4  and 
3.2  is  utilized  to  determine  the  minimal  required  silicon 
block/cantilever  width,  which  gives  rise  to  nano-tube  buck¬ 
ling.  This  width  was  found  to  be  14  silicon  lattice  parame¬ 
ters,  and  is  in  very  good  agreement  with  its  counterpart  (13 
silicon  lattice  parameters)  determined  using  molecular  me¬ 
chanics  simulations.  The  details  of  the  molecular  modeling 
simulations  can  be  found  elsewhere  [10]. 

To  further  validate  the  finite  element-based  approach  de¬ 
veloped  in  the  present  work,  a  comparison  is  made  in  Fig.  5 
between  the  total  strain  energies  (the  strain  energies  stored  in 
the  silicon  cantilever  nano-beam  and  the  SWCNT)  as  a  func¬ 
tion  of  the  cantilever  curvature  obtained  using  the  present 
approach  and  molecular  modeling.  It  is  seen  that  the  agree¬ 
ment  between  the  two  sets  of  results  is  good  both  relative  to 
the  curvature  dependence  of  the  strain  energy  and  relative  to 
the  value  of  critical  curvature,  at  which  nano-tube  buckling 
occurs. 

The  results  presented  in  the  present  section  are  encour¬ 
aging  and  suggests  that  the  present  finite  element-based 
approach  is  a  valuable  alternative  for  modeling  the  per¬ 
formance  of  functionalised  beam/carbon-nano-tube-based 
nano-devices. 
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3.2.  Optimization  of  the  nano-valve  design 

In  this  section,  the  finite  element-based  method  developed 
in  Section  2.6  is  used  to  identify  the  optimal  geometry  of  the 
nano-valve  under  consideration. 

The  first  step  toward  designing  the  nano-valve  is  to  select 
the  appropriate  SWCNT.  As  explained  earlier,  it  is  desirable 
to  have  a  SWCNT  with  a  sufficiently  large  diameter  to  enable 
an  easy  passage  of  the  molecules,  when  the  nano-valve  is  in 
the  open  position.  However,  previous  molecular  mechanics 
studies  e.g.  [10]  have  shown  that  SWCNTs  can  maintain  a 
stable  circular  cross-section  only  up  to  a  radius  of  approxi¬ 
mately  3  nm.  In  the  (n,  n)  arm-chair  family  of  SWCNTs,  the 
nano-tube  with  the  largest  radius  capable  of  retaining  a  sta¬ 
ble  circular  cross-section  is  found  to  be  a  (15,  15)  SWCNT. 
In  the  arm-chair  SWCNTs  with  a  larger  radius,  the  circular 
cross-section  subjected  to  compression  (in  a  direction  normal 
to  the  nano-tube  axis)  may  collapse  to  a  shape  consisting  of 
a  central  region,  in  which  the  opposite  walls  of  the  nano-tube 
are  nearly  parallel  to  each  other  and  at  a  distance  compara¬ 
ble  to  the  basal-plane  spacing  in  graphite.  The  side  portions 
of  the  cross-section  consist  of  nearly  circular  regions  with 
a  radius  of  approximately  0.75  nm,  which  appears  indepen¬ 
dent  of  the  value  of  the  chiral  number  n  >  15  in  the  arm-chair 
SWCNTs.  In  the  arm-chair  family  of  SWCNTs  between  (16, 
16)  and  (45,  45),  either  the  circular  or  the  collapsed  cross- 
sections  may  be  stable  depending  on  the  magnitude  of  the 
compression  applied  to  the  SWCNT  with  an  initial  circular 
cross-section.  At  low-levels  of  the  compression,  the  circu¬ 
lar  cross-section  is  stable  while  at  large-levels  of  the  com¬ 
pression,  the  collapsed  cross-section  is  more  stable.  In  the 
arm-chair  family  of  SWCNTs,  nano-tubes  with  the  chirality 
index  n  >  45,  only  the  collapsed  cross-section  becomes  stable. 
In  other  words,  no  energy  barrier  exist  in  the  configurational 
space  exist  between  the  circular  and  collapsed  cross-sections 
of  the  SWCNTs. 


Fig.  5.  Variation  of  the  total  strain  energy  with  the  silicon  nano-beam  cur¬ 
vature. 


A  stability  analysis  for  the  circular  cross-section  of  the 
SWCNTs  is  carried  out  in  this  work  using  the  present  finite 
element  formulation.  The  results  obtained  show  that,  for  the 
arm-chair  family  of  SWCNTs,  the  nano-tubes  with  the  chi¬ 
rality  index  n  <  17  have  a  stable  cross-sectional  shape,  while 
those  with  n  >  43  are  stable  in  the  collapsed  cross-sectional 
shape.  For  the  SWCNTs  with  an  intermediate  value  of  n ,  ei¬ 
ther  of  the  two  cross-sections  can  be  stable,  depending  on 
the  type  and  magnitude  of  the  perturbation  applied  to  the 
nano-tube  with  a  circular  cross-section.  The  observed  good 
agreement  between  the  results  of  the  atomic-level  and  the 
continuum-level  analyses  of  stability  of  the  circular  cross- 
section  in  the  SWCNTs  suggests  the  finite  element  formula¬ 
tion  proposed  in  the  present  work  can  realistically  represent 
deformation  behavior  of  the  SWCNTs. 

A  loss  of  the  circular  cross-section  is  not  acceptable  in  the 
present  case,  since  it  prevents  the  passage  of  large  molecules 
and  interferes  with  the  normal  “ on/off  ’  switching  of  the  nano¬ 
valve.  Consequently,  a  (17,  17)  SWCNT  with  a  radius  of 
1.15  nm  is  selected.  Such  SWCNT  offers  a  compromise  be¬ 
tween  a  desirable  large  cross-section  area,  required  stabil¬ 
ity  of  the  SWCNT  circular  cross-section  and  the  SWCNT’s 
ability  to  recover  a  large  (post-buckling)  localized  strain. 
The  length  of  the  (17,17)  SWCNT  is  next  arbitrarily  set  to 
~50nm.  To  match  the  length  of  the  SWCNT,  the  lengths  of 
the  silicon  block  and  the  cantilever  nano-beam  are  each  next 
set  to  ~20  nm.  Furthermore,  a  sufficiently  large  width  for  the 
silicon  block/cantilever  nano-beam  of  ~  108.6  nm  (20  silicon 
lattice  parameters)  is  selected,  so  that,  when  the  charge  on 
each  functional  group  of  the  monolayer  molecules  is  chosen 
arbitrarily  as  0.3e  (e  is  the  unit  electrical  charge),  SWCNT 
buckling  can  take  place  during  cantilever  actuation. 

Using  the  aforementioned  values  for  the  nano-valve  com¬ 
ponents,  a  finite  element  simulation  run  of  the  nano-valve 
closing  operation  due  to  cantilever  actuation  is  carried  out 
up  to  the  point  of  nano-tube  buckling.  The  morphology  of 
the  nano-valve  at  the  onset  of  SWCNT  buckling  is  shown  in 
Fig.  5.  Following  the  procedure  presented  in  Section  2.4,  the 
distribution  of  the  initial  required  nodal  forces  along  the  can¬ 
tilever  length  is  next  determined.  This  distribution  is  shown 
in  Fig.  4,  the  solid  curve  labeled  “minimal  required  initial 
force”.  Next,  a  set  of  dashed  curves  displayed  in  Fig.  4  is 
generated,  each  corresponding  to  the  initial  distribution  of  the 
nodal  forces  along  the  cantilever  length  in  a  nano-valve  of  a 
given  width.  To  meet  the  requirements  that  the  initial  nodal 
force  is  higher  than  the  required  minimal  force  for  each  node 
and  that  the  nano-valve  width  is  minimal,  the  lowest  dashed 
curve,  which  is  above  the  solid  curve  (the  one  labeled  “mini¬ 
mal  required  initial  force”)  is  selected  to  identify  the  minimal 
required  nano-valve  width  (~7.06nm=  13  silicon  lattice  pa¬ 
rameters). 

Once  the  minimal  required  width  of  the  nano-valve  is  de¬ 
termined,  the  finite  element  method  is  applied  again  for  the 
redesigned  nano-valve  to  verify  that  the  device  is  indeed  func¬ 
tional.  Furthermore,  a  device  with  a  width  of  6.52  nm  (i.e. 
one  silicon  lattice  parameter  smaller  width  than  the  minimal 


60 


M.  Grujicic  et  al.  /  Materials  Science  and  Engineering  B  117  (2005)  53-61 


Fig.  6.  An  example  of  the  nano- valve  configuration  at  the  onset  of  SWCNT 
buckling. 

required  width)  is  also  analyzed  using  the  finite  element 
method  to  verify  that  such  a  nano-device  would  not  func¬ 
tion,  since  the  reduction  of  the  electrostatic  energy  during 
cantilever  nano-beam  bending  is  not  sufficient  to  compen¬ 
sate  for  the  elastic  strain  energy  build-up  in  the  nano-beam 
and  in  the  SWCNT. 

The  geometrically  optimized  nano-valve  described  above 
contains  94,710  atoms  including  518  functionalized  groups 
in  the  top  face  monolayer.  Molecular  modeling,  for  the  rea¬ 
sons  presented  in  Section  1,  is  not  suitable  for  the  analysis  of 
a  device  with  this  many  of  atoms.  Solares  et  al.  [3]  proposed  a 
hybrid  method  within  which  the  individual  components  of  the 
nano-device  are  modeled  using  molecular  modeling,  while 


the  analysis  of  the  device  as  a  whole  is  carried  out  using  a  clas¬ 
sical  mechanical  design  approach.  We  applied  the  approach 
of  Solares  et  al.  [3]  to  the  geometrically  optimized  nano-valve 
design.  We  found  out  that  the  present  design  would  be  also 
functional  according  to  the  method  of  Solares  et  al.  [3].  How¬ 
ever,  the  method  of  Solares  et  al.  [3]  also  predicts  that  as  a 
device  with  a  width  smaller  by  two  silicon  lattice  parameters 
than  the  minimal  required  width  would  also  be  functional  in 
sharp  contrast  with  the  predictions  of  present  finite  element 
approach.  Since  both  the  hybrid  approach  of  Solares  et  al.  [3] 
and  the  present  finite  element  approach  utilize  atomic-level 
derived  properties  of  the  nano-valve  components,  the  dis¬ 
crepancy  between  the  two  approaches  can  be  explained  by 
recognizing  that  within  the  present  finite  element  approach, 
the  interactions  between  nano-valve  components  are  taken 
into  account,  while  such  interactions  are  not  considered  in 
the  approach  of  Solares  et  al.  [3].  An  example  of  the  ef¬ 
fect  of  the  interaction  between  the  cantilever  beam  and  the 
SWCNT  fluid-flow  conduit  is  seen  in  Fig.  6,  where  the  nano¬ 
tube  undergoes  a  localized  dimpling  deformation  at  the  point 
of  contact  with  the  cantilever  tip.  Such  deformation  is  asso¬ 
ciated  with  an  additional  strain  energy,  which  is  accounted 
for  in  the  present  work  but  not  in  the  approach  of  Solares  et 
al.  [3], 

3.3.  Design  considerations  for  an  “in-line”  nano-valve 

The  nano-valve  configuration  shown  in  Fig.  1  is  gener¬ 
ally  referred  to  as  “free-end”  design,  in  which  the  nano-tube 
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Fig.  7.  A  comparison  of  the  cross-sections  of  infinitely  long  (20,20)  arm-chair  SWCNTs  interacting  with  a  (100)  silicon  substrate  obtained  using:  (a)  an 
atomic-level  energy  minimization  analysis  and  (b)  a  finite  element  perturbation  analysis. 
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fluid-flow  conduit  can  freely  bend  under  the  effect  of  the 
cantilever  beam  until  the  point  of  buckling  is  attained.  There 
is  an  alternative  design  of  the  nano-valve  generally  referred 
to  as  an  “in-line”  design  [3].  In  this  design  of  the  nano¬ 
valve,  a  SWCNT  fluid-flow  conduit  is  positioned  on  the  top 
of  a  rigid  silicon  single-crystal  substitute  (beam).  Conse¬ 
quently,  bending  of  the  silicon  cantilever  during  the  clos¬ 
ing  operation  of  the  nano-valve  causes  local  crimping  of 
the  nano-tube,  which  stops  the  flow  of  the  fluid.  Since 
the  SWCNT  lies  on  the  silicon  substitute,  its  atoms  inter¬ 
act  with  those  of  the  substitute,  causing  a  change  in  the 
SWCNT  cross-section.  A  severe  departure  of  the  SWCNT 
cross-section  from  a  circular  shape,  when  the  nano-valve 
is  in  the  open  position  is  generally  not  desirable,  since  it 
may  interfere  with  the  flow  of  the  fluid.  In  this  section, 
we  conducted  a  test  of  the  present  finite  element  method 
to  determine  if  this  method  can  correctly  account  for  the 
SWCNT/substrate  interaction  induced  loss  of  the  SWCNT 
circular  cross-section. 

A  comparison  of  the  atomic-level  and  the  finite  element 
method  results  of  the  optimal  cross-section  shape  of  an  in¬ 
finitely  long  (20,  20)  SWCNT  in  contact  with  a  silicon  sub¬ 
stitute  is  given  in  Fig.  7(a)  and  (b).  A  fairly  good  agreement 
between  the  two  sets  of  results  shown  in  Fig.  7(a)  and  (b) 
suggest  that  the  present  finite  element  approach  can  be  used 
for  the  design  of  an  in-line  type  of  the  nano-valve.  A  detailed 
design  analysis  of  an  in-line  nano-valve  will  be  carried  out 
in  our  future  work. 

4.  Conclusions 

Based  on  the  results  obtained  in  the  present  work,  the 
following  main  conclusions  can  be  drawn: 

1 .  A  finite  element  method-based  procedure  for  the  analysis 

of  nano-size  devices,  in  which  doubly  curved  shell  el¬ 
ements  are  used  to  represent  the  SWCNT  wall  yield  the 


results,  which  is  quite  comparable  to  their  molecular  mod¬ 
eling  counterparts. 

2.  The  use  of  a  continuum-based  procedure,  such  as  the  one 
developed  in  the  present  work,  enables  the  analysis  of 
nano-size  devices,  for  which  a  purely  molecular  analysis 
would  be  either  impractical  or  impossible. 
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